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ABSTRACT 


Mean  turbulent  velocity  distribution  in  a  conduit 
of  rectangular  cross  section  was  determined  with  a  pitot 
tube.  The  aspect  ratio  of  the  test  section  was  8:1,  the 
width  8  inches  and  the  height  1  inch.  The  pitot  tube  used 
was  constructed  from  hypodermic  syringe  needles  and  the 
tip  outside  diameter  was  0.016  inches.  Tlie  inside  diameter 
of  the  pitot  tube  was  0.008  inches  and  the  static  tap  was 
0.0625  inches  in  diameter.  The  experiments  were  carried 
out  with  water  as  the  fluid. 

Velocity  distributions  were  determined  24  feet  or 
162  equivalent  diameters  from  the  upstream  end  for  water 
flow  rates  corresponding  to  Reynolds  numbers  of  5,070,  11,700, 
13,650  and  20,500.  An  attempt  to  correlate  the  observed 
data  with  u'*’  -  y"**  parameters  was  made. 

Secondary  flow  was  found  to  exist  in  the  range  of 
Reynolds  numbers  investigated.  The  relative  magnitude  of 
the  secondary  flow  as  compared  to  the  main  flow  seems  to 
decrease  with  increasing  Reynolds  number. 

Tlie  effect  of  secondary  flow  was  felt  even  at  the 
centre-line  velocity  distribution  in  the  channel  used  for 


the  investigation. 
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A  computer  programme,  based  on  the  theoretical 
method  proposed  by  Deissler  and  Taylor,  was  written  to 
calculate  the  velocity  distribution.  Since  this  method 
does  not  consider  the  existence  of  secondary  flow,  it  may 
be  used  to  provide  the  undisturbed  velocity  distribution 
in  the  conduit.  This  may  be  then  used  as  the  basis  for 
characterizing  the  secondary  flow. 


ACKNOWLEDGEMENTS 


The  author  wishes  to  express  his  sincere  gratitude 
to  Dr.  L.U.  Lilleleht  for  his  guidance  and  helpful  criticism 
during  this  investigation. 

Thanks  are  due  to  the  staff  of  the  department  work¬ 
shop  for  their  help  and  co-operation  in  constructing  the 
equipment. 

The  author  also  thanks  Mrs,  Georgina  Smith  for 
typing  the  thesis. 

The  financial  support  of  the  National  Research  Council 
of  Canada  is  gratefully  acknowledged. 


TABLE  OF  CONTENTS 


Page 

List  of  Tables  i 

List  of  Figures  ii 

INTRODUCTION  1 

LITERATURE  SURVEY  5 

Velocity  Measurement  Instruments  and  Techniques  5 

Velocity  Distribution  in  Rectangular  Conduits  15 

Velocity  Distribution  Correlations  16 

Entrance  Length  for  the  Development  of  Turbu¬ 
lent  Flov/  19 

EXPERIMENTAL  EQUIPMENT  20 

Rectangular  Conduit  20 

Test  Section  20 

Traversing  Mechanism  21 

Pitot  Tube  22 

Auxiliary  Equipment  22 

Manometer  24 

EXPERIMENTAL  PROCEDURE  25 


Page 


RESULTS  AND  DISCUSSION  27 

Preliminary  Experiments  with  Tube  27 

Computer  Programme  30 

Channel  Experiments  34 

CONCLUSIONS  AND  RECOMMENDATIONS  41 

BIBLIOGRAPHY  44 

APPENDICES  47 

I  Calibration  of  the  Pitot  Tube  48 

II  Computer  Programme  49 

III  Data  \.\rith  Tube  -  Preliminary  Experiments  52 

IV  Data  with  Channel  56 

V  u'*’  -  y"*"  Calculations  72 


NOMENCLATURE 


90 


1. 


LIST  OF  TABLES 

Pages 


Tables 

1 

to 

6 

Data  with  Tube 

53  - 

55 

Tables 

7 

to 

10 

Data  with  Channel 

56  - 

71 

Table 

11 

Centre  Line  u'*'  -  y'*’ 

Data 

72 

Table 

12 

u+  -  y'*"  Data  1  Inch 

from  Centre 

73 

1.1 


LIST  OF  FIGURES 

Page 

Figure  1  Schematic  Diagram  of  the  Experimental 

Equipment  74 

Figure  2  The  Traversing  Mechanism  75 

Figure  3  u**"  -  y'*’  Plot  of  the  Data  on  Tube  7G 

Figure  4  Velocity  Distribution  -  Method  of  Deissler 

and  Taylor,  Graphical  Computation  77 

Figures  5  Velocity  Distribution  along  Vertical 

to  10  Traverses  at  a  Reynolds  Number  of  13,650  73  -  83 

Figure  11  Constant  Velocity  Lines  at  a  Reynolds 

Number  of  13,650  84 

Figure  12  Prandtl's  Picture  of  Secondary  Flow 

and  Velocity  Distribution  along  a  Vertical 
Traverse  85 

Figure  13  u"’"  -  y**"  Plot  of  Centre  Line  Traverses 

Data  86 

Figure  14  u"^  -  y"*"  Plot  of  Data  on  Traverses  1  Inch 

from  the  Centre  87 

Figure  15  Approximate  Shear  Stress  at  the  Wall  88 

Figure  16  Method  of  Deissler  and  Taylor 


89 


INTRODUCTION 


In  the  transportation  of  viscous  fluids,  such  as  for 
example  crude  petroleum,  the  high  pressure  drop  in  the  flow 
line  is  the  main  factor  contributing  to  the  high  power 
requirement  and  to  the  consequent  cost  of  transportation. 

When  a  fluid  is  flowing  through  a  conduit  with  turbulent 
motion,  a  major  fraction  of  the  frictional  drag  is  concen¬ 
trated  in  a  very  thin  layer  of  the  fluid  close  to  the  conduit. 
If,  then,  a  film  of  a  less  viscous  fluid  is  substituted  in 
this  region,  the  frictional  drag  will  be  reduced.  Substitution 
of  this  resistance-film  of  petroleum  by  a  film  of  water  is 
certain  to  lessen  considerably  the  power  requirements  for 
transportation. 

Although  concentric  flow  of  the  two  immiscible  fluids 
is  more  advantageous  for  the  reduction  of  the  pressure 
gradient,  it  raises  problems  in  the  generation  and  maintenance 
of  such  a  configuration.  Even  a  slight  density  difference 
between  the  two  fluids  is  sufficient  to  cause  stratification 
unless  the  conduit  is  preferentially  wet  by  the  less  viscous 
phase.  This  makes  the  possibility  of  stratified  flow  more 
common  and  an  important  field  of  investigation. 
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In  an  earlier  project  at  the  University  of  Alberta 
Charles  '  ’  has  investigated  the  co-current  stratified  flow 
of  two  immiscible  liquids  in  a  rectangular  conduit.  In 
this  investigation  the  laminar- laminar  and  the  laminar- 
turbulent  case  were  studied.  Quantitative  expressions  were 
developed  for  the  velocity  distribution  and  the  volumetric 
flow  rate  for  the  laminar- laminar  case  only.  The  lack  of 
turbulent  velocity  distribution  data  in  rectangular  channels 
did  not  permit  the  development  of  similar  expressions  for 
the  laminar- turbulent  or  the  turbulent- turbulent  case. 

The  present  investigation  is  a  step  towards  the  quan¬ 
titative  study  of  stratified  flow  with  at  least  one  of  the 
phases  in  turbulent  motion.  This  includes  the  two  cases  of 
laminar- turbulent  and  turbulent- turbulent  flow,  in  order 
to  understand  the  phenomena  of  two- phase  flow  better,  it  is 
necessary  to  understand  the  flow  of  a  single  phase  in  turbulent 
motion  in  rectangular  channels. 

In  the  literature  on  the  mechanics  of  fluids,  most 
of  the  general  expressions  describing  the  velocity  distri¬ 
bution  have  been  obtained  for  the  geometrically  simpler 
cases  of  either  (i)  infinite  parallel  plates  or  (ii)  conduits 
of  circular  cross  sections.  Both  these  cases  are  amenable 
to  a  mathematical  approach  and  various  theories,  such  as 
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Prandtl's  mixing  length  theory,  von  Karman's  similarity 
theory,  Taylor's  statistical  approach  have  been  proposed. 

These  have  resulted  in  certain  theoretical  correlations  like 
the  logarithmic  velocity  distribution  law  and  semi-empirical 
correlations  like  the  expressions  of  Nikuradse  or  Deissler 
for  the  universal  velocity  distribution  law.  However,  the 
model  of  infinite  parallel  plates  is  true  only  in  theory. 

In  practice  it  is  approximated  by  flow  in  rectangular  channels 
with  high  aspect  ratios  so  that  the  situation  near  the  cen¬ 
tral  region  can  be  considered  to  be  the  same  as  that  between 
infinite  parallel  plates.  For  the  flow  field  across  the 
entire  cross  section  of  a  rectangular  conduit,  however,  the 
effect  of  the  side  walls  cannot  be  ignored. 

The  aim  of  this  investigation  was  primarily  to 
determine  the  velocity  distribution  in  turbulent  flow  for  a 

(23) 

conduit  of  a  rectangular  cross  section.  Earlier  data  of  Nikuradse'  ' 
showed  a  tendency  for  the  constant  velocity  lines  to  shift 
towards  the  corners.  This  had  been  explained  by  Prandtl  (26) 
to  be  caused  by,  what  he  described  as,  "secondary  flow". 

It  was  also  the  aim  of  this  investigation  to  determine  the 
extent  of,  and  if  possible  to  characterize,  the  secondary  flow. 

The  experiments  were  carried  out  in  a  rectangular 
conduit  with  an  aspect  ratio  of  8:1,  The  velocity  distri- 
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bution  was  determined  by  a  pitot  tube  using  a  manometer  with 
a  small  density  difference.  All  the  runs  were  made  with 
water,  and  the  results  can  be  extended  to  any  other  incom¬ 
pressible  fluid  using  the  concept  of  dynamic  similarity. 
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LITERATURE  SURVEY 

In  this  section  some  of  the  point  velocity  measurement 
devices  will  be  reviewed,  and  their  merits  and  demerits  with 
regard  to  the  present  investigation  will  be  considered. 

A  brief  review  of  some  of  the  velocity  distribution  laws  and 
correlations  will  follow  and  finally  the  minimum  entrance 
length  requirement  will  be  discussed. 

Velocity  Measurement  Instruments  and  Techniques; 

For  the  point  measurements  of  velocities  in  a  conduit  the 
following  are  some  of  the  instruments  and  techniques  available: 

(a)  Flow  visualization 

(b)  The  hot-wire  anemometer 

(c)  The  hot  film  anemometer 

(d)  Electromagnetic  velometry 

(e)  The  pitot-static  tube 

(f)  The  pitot  tube  with  a  separate  static  tube. 

(a)  One  of  the  oldest  techniques  of  measurements  is  based 
on  making  the  flow  pattern  visible  or  perceptible  to  a  suitable 
detecting  element,  like  a  photographic  plate.  This  is  the 
flow  visualization  technique.  In  his  classic  experiments, 
Osborne  Reynolds  used  dyes  for  visually  observing  the  flow. 
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With  gas  flow  wood  smoke,  hexite,  vapours  of  ammonia,  iodine 
etc.  have  been  used.  Instead  of  using  such  continuous  materials 
even  discrete  particles  can  be  introduced  in  the  flow  so 
long  as  the  size  of  the  particles  is  smaller  than  the  micro 
scale  of  turbulence  and  the  difference  between  the  density 
of  the  particles  and  the  fluid  density  is  not  great.  For 
liquids  aluminum  dust  and  emulsions  of  other  liquids  have 
been  used  besides  the  dyes.  With  air  tiny  soap  bubbles  can 
be  used.  The  path  of  the  discrete  particles  is  then  followed 
by  cinematography  or  stroboscopic  flash  photography.  Charles 
in  his  investigation  of  two  phase  flow  has  used  aluminum 
powder  with  stroboscopic  flash  photography  for  the  study  of 
the  interfacial  waves. 

(b)  The  hot-wire  anemometer  has  been  used  quite  extensively 
for  work  in  the  turbulent  region  (H' 27)  ^  It  uses  the  principle 
of  change  of  resistance  of  a  conductor  with  temperature  for 
the  measurement  of  velocities.  As  the  fluid  flows  faster, 
rate  of  heat  transfer  increases  and  the  change  in  resistance 
of  the  wire  due  to  the  drop  in  temperature  can  be  measured 
very  accurately  on  a  Wheatstone  bridge  circuit.  The 
anemometer  can  be  operated  at  (i)  constant  current  (ii)  con¬ 
stant  resistance  or  (iii)  constant  voltage.  Hinze^^^^ 
mentions  that  the  constant  current  method  has  been  developed 
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into  a  nearly  foolproof  tool,  whereas  the  electronic  system 
of  the  constant  resistance  method  is  not  as  reliable. 

Prandtl  and  Tietjens^^^^  mention  the  constant  voltage  opera¬ 
tion  and  indicate  its  suitability  to  be  only  for  very  small 
velocities.  A  detailed  description  of  the  constant  current 
and  constant  resistance  operations  can  be  found  in  Hinze . 

The  hot-wire  anemometer  is  a  very  sensitive  instrument 
capable  of  measuring  the  turbulent  fluctuations  in  a  flow. 

With  suitable  electronic  equipment  to  compensate  for  the 
non-linear  characteristics  of  the  circuit  and  the  wire  lag, 
it  is  possible  to  obtain  a  true  picture  of  the  turbulent 
fluctuations  on  a  cathode-ray  oscilloscope.  Various  arrange¬ 
ments  and  methods  have  been  developed  to  measure  turbulence 
characteristics,  and  by  suitable  manipulations  it  is  now 
possible  to  measure  the  turbulence  intensity  in  any  spatial 
direction,  the  turbulent  shear  component,  or  Reynolds'  shear 
stresses  and  u'w*",  the  Eulerian  correlation  coefficient, 

the  scale  of  turbulence  and  the  energy  spectrum  of  turbulence. 

However,  the  size  of  the  wire,  which  is  usually  0.0001 
to  0.0003  in.  in  diameter  and  0,025  to  0,06  in.  in  length 
restricts  the  use  of  the  hot-wire  anemometer  to  the  flow  of 
clean  gases.  The  sensitivity  is  affected  by  deposits  of  dust 
which  change  the  heat  transfer  characteristics  of  the  wire. 
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Any  heavier  particles  of  dust  can  cause  damage  to  the  wire. 
Frequent  cleaning  of  the  wire  is  essential,  as  is  the  checking 
of  the  calibration  curve. 

Frenkiel^^^  has  discussed  the  effects  of  the  wire 
length  in  the  anemometer.  A  correction  for  hot  wire  readings 
in  a  boundary  layer  for  proximity  to  the  solid  boundary  has 
been  suggested  by  Livesey^^^^. 

Kovasznay has  developed  a  hot-wire  turbulence 
measuring  equipment  suitable  for  measuring  accurately  the 
turbulent  fluctuations  at  supersonic  velocities  with  gases. 

The  equipment  covers  a  frequency  range  from  2  to  about  70,000 
cycles  per  second  and  uses  the  constant  current  mode  of 
operation. 

The  hot-wire  anemometer  finds  very  little  application 
with  liquids.  A  low  operating  temperature  is  necessary  in 
the  case  of  liquids  because  of  the  possibility  of  vapouriza¬ 
tion  and  the  formation  of  scale.  Middlebrook  and  Piret^^^^ 
have  pointed  out  the  possibility  of  deposits  forming  even 
in  rather  pure  liquids.  Hence  a  sturdier  construction  with 
thicker  and  longer  wires  is  necessary,  which  affects  the 
sensitivity  and  the  time  constant  of  the  wire  adversely. 

Longer  wires  are  undesirable  also  because  of  their  greater 
spatial  averaging  effect.  Hot  wire  anemometers  have  been 
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used  with  liquids  for  measurements  of  mean  velocities  in 
boundary  layers,  for  the  study  of  transition  phenomena  or 
turbulence  of  rather  low  frequencies. 

A  summary  report  on  the  development  of  hot-wire  tur¬ 
bulence  sensing  element  for  use  in  water  has  been  published 
by  Stevens,  Borden  and  Strauser .  The  report  covers  the 
work  done  between  1946  and  1956.  The  use  of  alternating 
current  has  been  recommended  for  work  with  water  and  the 
difficulties  regarding  cleaning  and  calibration  when  working 
with  ordinary  water  have  been  noted.  Vernotte  criticises 
the  use  of  the  law  of  heat  transfer  used  in  the  derivation 
of  the  correlation  between  the  temperature  and  the  current 
for  the  hot-wire  anemometer ,  The  law  used  is  valid  only 
under  steady  conditions,  whereas  in  turbulence  measurements 
the  wire  is  exposed  to  a  fluid  flow  of  unsteady  character. 

(c)  Ling  and  Hubbard have  developed  a  hot-film  anemometer 
for  use  with  liquids.  The  delicate  nature  of  the  hot-wire 
anemometer  is  due  to  the  thin  wire  stretched  across  the  two 
ends  of  a  prong.  In  the  hot  film  anemometer  the  wire  is 
replaced  by  a  film  of  platinum  on  a  chisel-shaped  glass  probe 
tip.  The  leads  from  the  film  are  taken  out  from  within  the 
probe.  The  operating  principle  is  the  same  as  that  for  the 
hot-wire  anemometer.  Since  the  mechanical  strength  of  such 
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a  film  is  much  higher  than  that  of  the  hot-wire  element,  it 
is  more  suitable  for  use  with  liquids.  However,  the 
restrictions  regarding  the  temperature  of  operation,  possibility 
of  formation  of  a  dirt  scale  and  the  objection  to  the  heat 
transfer  law  employed,  remain. 

(d)  Another  development  in  the  field  of  turbulent  measure¬ 
ments  is  the  technique  of  Electromagnetic  Velometry,  intro¬ 
duced  by  Kolin^^^.  This  method  employs  the  principle  of 
electromagnetic  induction.  According  to  Faraday's  law  of 
electromagnetic  induction  an  electromotive  force  is 
induced  in  a  conductor  moving  relative  to  a  magnetic  field. 

The  ionized  flowing  fluid  can  be  regarded  as  a  moving  conductor, 
in  which  case  the  induced  electric  field  strength  is  equal 
to  the  vector  product  of  the  magnetic  field  strength  and  the 
velocity.  Thus  if  U  is  the  velocity  component  perpendicular 
to  an  electromagnetic  field  of  strength  H,  the  strength  V 
of  the  induced  electric  field  is 

V  =  — HU 

c 

where  |i  is  the  relative  magnetic  permeability  of  the  liquid 
and  c  is  the  velocity  of  light. 

The  probe  is  in  the  form  of  two,  closely  spaced,  fine 
wires.  If  the  distance  between  the  wires  is  s,  the  potential 
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difference  E  will  be 

E  =  Vs  =  HUs 

c 

This  impulse  can  be  amplified  and  the  root  mean  square  measured 
by  a  thermal  milliammeter.  Hinze  points  out  the  importance 
of  the  linearity  of  the  voltage  with  the  velocity  component 
perpendicular  to  the  connecting  line  through  the  electrodes 
and  to  the  direction  of  the  electromagnetic  field. 

Grossman  et  have  shown  that  there  is  no  thermal 

lag  in  this  instrument  and  unlike  the  hot  wire  anemometer  its 
use  is  not  restricted  to  gas  flows.  Since  the  transient 
potentials  are  of  the  order  of  lO"*^  to  10“^  volts,  fairly 
elaborate  precautions  must  be  taken  against  any  spurious  potentials 
such  as  for  example  background  noise.  Use  of  non-conducting 
pipes  is  definitely  recommended.  Near  the  walls  a  large  velocity 
gradient  may  cause  "short-circuiting"  currents  which  in  turn 
may  modify  the  original  flov;.  However,  for  an  essentially 
flat  velocity  profile,  this  effect  is  negligible.  Kolin  and 
Reiche  have  analyzed  the  distribution  of  flow  through  a  homo¬ 

geneous  magnetic  field  perpendicular  to  the  floVv^  and  shown 
that  it  is  possible  to  eliminate  the  difficulties  arising 
from  induced  currents  by  measuring  the  potential  gradient 
along  a  diameter  bisecting  the  angle  between  the  H  axis 
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and  the  Y  axis  ^^rhere  the  latter  is  perpendicular  to  both  the 
H  and  the  U  axes. 

(e)  The  pitot- static  tube  was  not  used  in  the  present  investi¬ 
gation  considering  the  difficulty  in  the  manufacture  of  very 
small  diameter  tubes.  To  enter  the  buffer  zone  or  the  laminar 
sublayer  in  a  channel  of  1  inch  height,  the  pitot- static 
tube  must  be  of  a  diameter  less  than  0,016  inc.  This  is  easily 
calculated  assuming 


for  the  laminar  sublayer  and  using  the  concept  of  equivalent 
diameter  for  the  channel  with  turbulent  flov/  at  a  Reynolds 
number  of  4000.  Except  for  the  relative  insensitivity  of  the 
pitot-static  tube  to  the  angle  of  yaw^^^  ,  it  has  no  distinct 
advantage  over  the  simple  pitot  tube  with  a  static  opening 
built  in  the  wall. 

(f)  Among  the  instruments  and  techniques  of  point  velocity 
measurements  mentioned  above,  all  ?out  the  pitot- tube  and  the 
pi tot- static  tube,  are  mainly  used  for  measurements  of  turbulence 
quantities  such  as  the  fluctuating  velocity  components,  the 
Reynolds'  stresses,  the  correlation  coefficient,  etc.  For  the 
purpose  of  determining  the  mean  velocity  at  any  point,  the 
sensitivity  of  the  anemometers  or  the  electromagnetic  probe 
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will  be  a  drawback  rather  than  an  asset.  A  careful  averaging 
of  the  fluctuating  velocities  will  be  necessary. 

In  this  respect  then,  the  simple  pitot  tube  with  a 
static  hole  in  the  wall  seems  to  be  the  most  advantageous. 

The  natural  inertia  of  the  pitot- tube  -  manometer  system  causes 
a  time  average  value  to  be  registered.  The  relative  insensitivity 
prevents  the  turbulent  fluctuations  from  being  indicated. 

Among  other  advantages  are  the  simple  construction,  necessity 
of  little  auxiliary  equipment  and  relatively  easy  measurements. 

The  calculations  to  obtain  the  velocity  distribution  are  also 
straight-forward, 

Folsom has  presented  an  excellent  review  of  the 
available  literature  on  the  Pitot  Tube.  Hinze has  dis¬ 
cussed  the  difficulties  arising  in  the  accurate  measurement 
of  the  mean  static  pressure  in  a  turbulent  flow.  The  assum¬ 
ption  of  a  constant  mean  static  pressure  over  a  cross  section 
in  a  straight  duct  has  been  found  to  be  satisfactory,  Knudsen 
and  Katz  have  quantitatively  discussed  the  disadvantage 
of  using  a  pitot  tube  for  measurements  in  the  laminar  sublayer. 

It  can  be  shown  that  the  pressure  measured  by  a  pitot  tube 
when  very  near  the  wall,  may  correspond  to  the  velocity  some 
distance  away  from  the  actual  centre  of  the  probe,  towards  the 
region  of  higher  velocity.  Such  a  difference  has  also  been 
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reported  by  Young  and  Maas^^^^.  The  possibility  of  circulation 
within  the  probe  tip  also  exists,  this  may  cause  still  higher 
pressure  readings  to  be  registered. 

Velocity  Distribution  in  Rectangular  Conduits; 

No  extensive  studies  have  been  conducted  on  the  velocity 
distribution  in  conduits  having  non-circular  cross  sections. 
Prandtl  and  Nikuradse  have  determined  point  velocities  in 
such  conduits  and  also  studied  the  flow  patterns  by  visual 
techniques  .  These  studies  revealed  the  existence  of  a 

flow  towards  the  corners  of  the  conduit  and  away  from  the 
sides,  superimposed  upon  the  main  longitudinal  motion.  The 
effect  of  this  "secondary  flow"  in  pushing  the  lines  of 
constant  velocity  towards  the  corner  and  away  from  the  v/all 
has  been  illustrated  in  a  later  section. 

Secondary  flow  associated  with  fully  developed  turbulent 
flow  are  reported  to  have  been  studied  by  Hoagland,  LcC. 

n  2) 

in  long,  straight  rectangular  ducts'  '  .  Recently  Shook,  C.A. 
and  Daniel,  S.M.  have  presented  a  paper  on  the  turbulent 

flow  velocity  distribution  on  the  axis  of  a  rectangular  channel. 
They  found  a  significant  deviation  from  Nikuradse ' s  results. 

The  investigation  was  carried  out  in  a  rectangular  channel 
0.97"  X  3.976"  in  cross  section  and  the  pitot  tube  V7as  0.058 


inches  in  diameter. 
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Velocity  Distribution  Correlations; 

One  of  the  earliest  velocity  distribution  correlations 
is  the  ^th  power  law.  On  the  basis  of  Nikuradse's  data  on 
velocity  profiles  in  smooth  pipes  it  can  be  noted  that  an 
empirical  relation  of  the  form 


u  _ 
U 


correlates  the  data  well.  The  value  of  the  exponent  n  is 

n  =  6  at  Re  =  4000;  it  increases  to  n  =  7  at  Re  ~  110  x  10^ 

and  to  n  =  10  at  Re  =  3240  x  10^  which  was  the  highest  Reynolds 

(29^ 

number  attained  in  the  investigation '  ‘  , 

Prandtl  discovered  a  relation  between  the  — th  power 

n 


lav7  and  Blasius'  law  of  friction 


f  =  0.0792  (Re) 

and  obtained  the  power  velocity  distribution  law 

JL.  =  3.74 

u* 

1/2 

where  u*  =  (tq/p)  is  the  friction  velocity. 

Since  the  time  Prandtl  put  forward  his  mixing  length 
theory  to  derive  the  logarithmic  velocity  distribution  law 
for  flov/  over  flat  plates,  a  nuraber  of  other  correlations 
have  appeared.  Most  of  these  investigators  seem  to  accept  the 
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model  of  turbulent  flow  consisting  of  a  central  turbulent  core 
and  a  laminar  sublayer  next  to  the  solid  surface  with  a 
transitional  or  buffer  zone  between  the  two.  The  laminar 
sublayer  can  be  represented  by 

u-"  =  y^<5 

Nikuradse  /  as  a  result  of  his  work  on  tubes,  gives 
the  following  relation 


y’^>30 

Deissler^^^  developed  a  single  equation  for  tubes 
relating  U''*  to  y**"  up  to  y"^  =  26,  For  the  turbulent  core 
this  relation  is 

U’’*  =  3.8  +  i —  In  y"^  y***  >  26 

0.360 

Skinner's  relation with  his  work  on  tubes  is 

=  6,0  +  — - —  In  y’’’ 

0.369 

Laufer^^^^  and  Wattendorf  made  measurements  with 
two  dimensional  flow  in  channels  and  came  up  with 


Laufer : 


5.5 


0 . 344 


1 


In  y+ 
In  y+ 


Wattendorf : 


4.0  + 


0.373 
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Ross  critically  reanalyzed  Nikuradse's  data  and  came 
up  with^^^^ 

=  5,6  +  5.6  log  y"**  yAw  ^0.13 

Rothfus  and  Monrad  found  that  versus 

y"*"  gives  a  better  graphical  correlation. 

Wang^^^^  developed  a  formula  for  velocity  distribution 
from  a  study  of  mixing  length  distribution  across  the  section 
of  channels.  This  formula,  though  rather  cumbersome  to  handle, 
checks  well  with  experiments  and  yields  the  same 
friction  formula  as  derived  by  von  Karman  and  Prandtl. 

Pai^^^^  has  derived  an  equation  for  the  velocity 
distribution  from  Reynolds'  momentum  equations. 

Gill  and  Scher ,  in  a  recent  paper,  have  derived  a 
continuous  velocity  distribution  based  on  an  arbitrary 
modification  of  Prandtl 's  mixing  length  expression.  This 
continuous  expression  converges  to  the  Hagen- Poiseuille 
equation  at  appropriate  Reynolds  numbers  without  subdividing 
the  flow  in  separate  regions. 

It  has  been  noticed  that  the  universal  velocity  distri¬ 
bution  equations  mentioned  earlier  show  a  systematic  deviation 
from  most  experimental  data.  Expressions  of  the  type 

=  A  +  B  log  y"*" 
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have  been  suggested  by  Millikan,  Reichardt  and  Hinze^^^^. 

The  correction  function  C (y/K)  is  graphically  represented  on 
the  basis  of  available  experimental  data.  Bogue  and  Metzner^^^ 
in  a  recent  paper,  have  suggested  a  correction  function  as  a 
function  of  y/R  and  the  Fanning  friction  factor. 

Entrance  Length  for  the  Development  of  Turbulent  Flow; 

No  reliable  relationship  is  available  for  predicting 
the  minimum  entrance  length  required  for  the  full  development 
of  turbulent  flow.  Experiments  by  Kirsten  and  Nikuradse  show 
that  the  transition  length  is  much  shorter  than  that  for  laminar 
flow  and  is  less  dependent  on  the  Reynolds  number 
Kirsten  mentions  the  transition  length  to  be  100  to  200  radii, 
Nikuradse  50  to  80  radii,  whereas  Latzko  gives  a  figure  of 
20  radii  for  a  Reynolds  number  of  20,000. 
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EXPERIMENTAL,  EQUIPMENT 

In  this  chapter  the  equipment  used  for  the  experimental 
work  will  be  described.  A  schematic  diagram  of  the  experi¬ 
mental  setup  appears  in  Figure  1. 

Rectangular  Conduit; 

The  conduit  has  a  rectangular  cross  section  with  an 
average  aspect  ratio  of  7.95:1  and  an  overall  length  of 
37.5  feet.  The  entry  section  has  been  specially  designed  for 
two  phase  work,  such  as  to  cause  as  little  mixing  at  the  inter¬ 
face  as  possible.  The  fluids  enter  this  section  against  a 
baffle  at  a  circular  cross  section  which  gradually  changes  over 
to  a  rectangular  cross  section  over  a  distance  of  two  feet. 
However,  for  the  present,  single  phase  experiments,  the  upper 
ent2ry-line  previously  used  for  oil  was  kept  closed, 

A  detailed  description  of  the  conduit  has  been  given 
by  Charles  . 

Test  Section; 

The  conduit  has  a  removable  test  section  24  feet  from 
the  upstream  end,  previously  used  for  photographic  purposes. 
This  section  was  modified  for  the  velocity  measurements,  A 
traversing  mechanism  was  installed  so  as  to  enable  velocity 
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determination  at  any  point  over  the  cross  section.  A  careful 
measurement  of  the  dimensions  of  the  test  section  revealed 
the  aspect  ratio  to  be  exactly  8:1  at  the  cross  section  over 
which  the  tip  of  the  pitot  tube  travelled. 

Traversing-  Mechanism: 

The  traversing  mechanism,  Figure  (2) ,  consisted  of  a 
horizontally  sliding  plate,  sliding  across  the  channel,  with 
a  vertical  micrometer  screw  mounted  on  it.  The  two  motions 
provided  the  probe,  mounted  at  the  end  of  the  micrometer 
screw,  full  access  to  any  point  across  the  cross  section;  the 
only  limitation  being  due  to  the  dimensions  of  the  probe. 

The  horizontal  plate  was  made  of  tool  steel  and  was 
1-1/2  inches  x  2  inches  x  1/2  inch.  It  was  dovetailed  to 
a  2  inch  width  at  the  bottom.  The  female  section  of  the  dove¬ 
tail  V7as  made  of  brass  and  had  bolts  on  one  side  which  could 
press  a  thin  brass  plate  on  the  sliding  plate.  The  female 
section  v/as  cemented  to  the  top  of  the  test  sec-fion  and  a  slot 
was  cut  through  the  female  section  of  the  dovetailed  sliding 
plate  and  also  through  the  top  of  the  test  section  so  as  to 
allow  the  probe  to  move.  The  gasket  was  a  1/8  inch  thick 
rubber  sheet  with  a  paste  of  graphite  in  grease  used  as  the 
lubricant.  The  gasket  with  this  lubricant  was  satisfactory 
for  preventing  leaks  and  yet  keeping  the  tool  steel  plate 
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capable  of  easy  movement. 

The  tool  steel  plate  carried  another  dovetail  slide 
with  its  housing.  A  micrometer  screw  was  fixed  at  the  top  and 
the  spindle  could  be  used  to  measure  the  movement  of  the  slide. 
This  slide  carried  the  pitot  tube  which  was  soldered  to  its 
head  with  the  top  taken  out  from  the  side  of  the  slide. 

The  pitot  tube  passed  through  a  rubber  sleeve  in  the 
tool  plate  slide,  through  the  slot  and  into  the  channel, 

Pitot  Tube; 

The  pitot  tube  was  constructed  from  hypodermic  syringe 
needles.  Four  needles  of  Stubbs  gages  13,  19,  24  and  27  were 
used  to  achieve  a  step-wise  reduction  in  the  diameter  of  the 
tube.  Needles  of  gages  18  and  19  were  found  suitable  in 
stiffness  for  the  support  and  the  needle  of  gage  27  was  used 
as  the  probe  tip.  This  provided  a  pitot  tube  with  an  outer 
diameter  of  0.016  inches  and  an  inside  diameter  of  0.008  inches. 
The  tip  of  the  pitot  tube  was  carefully  ground  to  a  smooth 
cone.  Measurements  within  0.008  inches  of  the  walls  v/ere 
possible . 

Auxiliary  Equipment; 

Water  was  stored  in  a  tank  52  inches  in  diameter,  75 
inches  in  height,  with  a  removable  cover.  It  was  made  of  14 
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gauge  galvanized  iron  sheet  and  had  a  capacity  of  about  650 
U.S.  gallons. 

The  pump  used  was  of  the  positive  displacement  gear 
type,  model  FL  manufactured  by  Jabsco  Pump  Co.,  with  a  capacity 
of  about  30  U.S.  gallons  per  minute  at  1750  R.P.M.  The  drive 
was  from  a  3  phase  two  horse  power  electric  motor  running  at 
1750  R.P.M. 

Due  to  the  continuous  recirculation  in  the  tank  through 
the  pump  and  the  channel,  the  temperature  of  the  water 
increased  gradually.  This  necessitated  purging  of  the  hot 
water  from  the  storage  tank.  This  could  be  done  conveniently 
through  two  bottom  outlets  in  the  tank.  Fresh  make  up  water 
could  be  admitted  through  three  equally  spaced  1/2  inch 
diameter  openings  on  the  periphery  of  the  tank  27  inches 
from  the  bottom.  Mixing  of  the  cold  water  with  the  water  in 
the  tank  was  thus  achieved  quite  rapidly. 

Flow  meters  used  were  of  the  rotameter  type.  Three 
rotameters  were  installed  with  a  view  to  covering  a  range  of 
flow  rates  from  0.3  U.S.  gallons  per  minute  to  30  U.S.  gallons 
per  minute.  In  the  present  investigation  for  the  Reynolds 
number  of  20,500  the  three  rotameters  were  used  in  parallel. 
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Manometer ; 

The  differential  pressure  measurements  were  made  with 
a  manometer  made  of  10  miTi  glass  tubing.  The  simple  two- limb 
manometer  was  considered  quite  adequate  for  the  purpose  pro¬ 
vided  a  suitable  manometric  liquid  was  used.  The  leads  from 
the  pitot  tube  and  the  static  tap  were  of  plastic  tubing. 

Instead  of  using  a  graduated  scale  behind  the  mano¬ 
metric  tubes,  for  a  high  accuracy  a  cathetometer  was  used. 

The  cathetometer  had  a  tripod  base  with  adjustable  screws  for 
levelling.  The  telescope  also  carried  a  leveller  and  an 
adjustable  screw  for  adjusting  it  parallel  to  the  horizontal. 

The  use  of  the  cathetometer  made  it  possible  to  read  a  level 
difference  with  an  accuracy  of  jro.005  cm. 

The  manometric  fluid  used  was  a  mixture  of  Benzene  and 
Carbon  Tetrachloride.  Both  these  liquids  have  a  reported 
solubility  of  0.08  parts  per  100  parts  at  20 ®C.  This  solubility 
was  considered  low  enough  to  neglect  a  change  in  the  density 
of  the  water  in  the  limbs  over  the  organic  liquid.  The 
specific  gravity  of  the  mixture  used  v/as  determined  experimen¬ 
tally  and  found  to  be  1.0895  at  77°F. 
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EXPERIMENTAL  PROCEDURE 

Water  was  introduced  into  the  tank  with  the  valves  at 
both  the  ends  of  the  channel  open  so  that  the  channel  was 
filled  gradually  as  the  level  of  the  water  rose  in  the  storage 
tank.  This  was  done  in  order  not  to  subject  the  channel  to 
the  stresses  caused  by  the  intermittent  surging  of  water  from 
the  tank  into  the  channel  when  the  main  return  valve  is  opened 
on  an  empty  channel. 

With  the  water  in  the  tank  at  the  desired  level,  the 
pump  was  started  and  the  flow  rate  adjusted  to  the  required 
value  by  a  suitable  manipulation  of  the  valves. 

The  pitot  tube  was  brought  to  the  desired  point  and  the 
system  was  allowed  to  reach  equilibrium.  Theoretically, 
equilibrium  would  be  reached  only  after  an  infinite  time  inter¬ 
val;  but  for  all  practical  purposes  equilibrium  was  assumed 
to  have  been  reached  if  the  manometer  did  not  show  any  detectable 
difference  between  two  observations  made  half  an  hour  apart. 
Depending  on  the  deviation  from  equilibrium  the  time  taken 
for  one  observation  ranged  from  1-1/2  to  4-1/2  hours. 

The  telescope  V7as  kept  as  close  to  the  manometer  as 
possible  to  minimize  the  error  due  to  any  possible  inclination 
from  the  horizontal.  It  was  also  placed  at  an  equal  distance 
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from  both  the  limbs  so  that  refocussing  was  never  necessary. 

After  each  vertical  traverse  the  slide  was  moved  to  the 
next  position.  At  least  six  traverses  were  made  on  one  half 
the  width  of  the  channel. 

Readings  were  attempted  for  each  traverse  at  a  minimum 
of  12  positions  with  4  spots  symmetrically  placed  about  the 
midpoint  of  the  vertical  traverse. 

The  temperature  of  the  water  was  not  allowed  to  fluctuate 
more  than  ^2®?  about  the  desired  mean.  From  time  to  time  the 
hot  water  was  purged  out  from  the  tank  and  fresh  water  was 
added  to  lower  the  temperature.  The  pump  was  kept  running 
during  this  operation  so  that  the  mixing  was  uniform  and  rapid. 
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RESULTS  AND  DISCUSSION 


Preliminary  Experiments  with  Tube; 

The  preliminary  experiments  with  the  vertical  tube 
served  the  two  purposes  of  (i)  establishing  the  adequacy  of 
the  pitot  tube  as  a  mean  velocity  measuring  device,  and  (ii) 
calibrating  the  pitot  tube. 

The  experimental  data  on  the  tube  are  given  in 
Appendix  III.  Complete  velocity  profiles  were  determined  for 
the  following  values  of  the  Reynolds  number  9150,  17,690  and 
22,260.  Observations  were  also  made  at  Reynolds  numbers  of 
2,700  and  34,300. 

The  complete  profiles  were  integrated  to  obtain  the 
value  of  the  total  flow  through  the  tube.  These  values  were 
found  to  be  within  tl.2%  of  the  values  observed  by  the  rota¬ 
meter  installed  in  the  flow  line.  Table  6  compares  the 
integrated  values  with  the  rotameter  values.  Within  the 
limits  of  experimental  accuracy,  a  pitot  tube  coefficient 
of  unity  could  be  used. 

Reproducibility  of  the  observed  data  was  checked  by 
repeating  some  observations.  In  all  cases  the  agreement 
between  the  two  observations  was  v/ithin  tl%- 
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A  plot  of  the  data  with  dimensionless  co-ordinates 
u'*’  and  is  given  in  Figure  3.  For  comparison  with  the 
observed  data,  correlations  given  by  Laufer,  Skinner, 

Deissler  and  Nikuradse  have  also  been  drawn.  As  mentioned  in 
an  earlier  section,  Laufer  and  Nikuradse  used  rectangular 
channels,  while  Skinner  and  Deissler  used  smooth  tubes  for 
their  experimental  work.  Obviously  there  is  a  notable  diver¬ 
gence  between  the  correlations  of  Laufer  and  Skinner  on  the 
one  hand  and  Deissler  and  Nikuradse  on  the  other,  A  bulk  of 
the  data  from  the  present  investigation  falls  between  the 
correlations  of  Skinner  and  Deissler.  The  data  at  a  Reynolds 
number  of  34,300  approach  Nikuradse ' s  curve  closely. 

In  spite  of  the  small  size  of  the  pitot  tube,  it  was 
not  found  possible  to  make  many  observations  in  the  buffer 
zone  and  the  laminar  sublayer.  For  example,  at  a  ^^eynolds 
number  of  17,690  assiiming  the  laminar  sublayer  to  be  for  values 
of  u"^  and  y"**  less  than  5,  the  thickness  of  the  sublayer  can 
be  calculated  to  be  0.00362  inches.  Assuming  the  limit  of 
the  buffer  zone  to  be  at 

=  12.85  y+  =  26 

the  upper  bound  of  the  buffer  zone  can  be  calculated  to  be  at 
a  distance  of  0.01883  inches  from  the  wall  for  the  3/4  inch 
inside  diameter  plastic  tube  used  for  the  experiments.  The 


fli  *pms£i 


IT'  ''•  ■  .a*}’ 

:itol  aikcfii^^  ' 


“  r  4^^ 


e/ict  no-  3&nni^3 


^  •  '-=-r  ,» 


.,.r  -  ..  ■  '  ' 

aWofiY^a  ^ 

>.  ^  V,;  ^ 

*f--  %|i  -v-mjL  ww  J 

"''■^‘riffi  ■’“^.j^;»°  i»i:‘»»->*-  »?gi^^ 
.-Sfni  .  «M.W#f3i«^'A„l£j&  l-Jgf 

W^. ...  ^ 


•;£^, 


‘  jt  wwmrn 

-  -  /SI.  ::.„‘^‘ 

Sfi  “(-f  oa  5a»6lii$iBa  ad-neo  sftos"  loSiiii^  ari-^ 


„.,  “  4.. 

fe ' 

ran 

,01  ixsv,  L 

0^  .ai„a™i,aq*»  arid-^oX 


'■'•!  , 


■  -w  r'i* '’^‘'*i)  I 


29 


pitot  tube  used  for  these  experiments  was  of  0.013  inch 
outside  diameter.  A  few  observations  made  in  the  buffer 
zone  were  found  to  deviate  considerably  from  Nikuradse ' s 
curve.  The  observations  tended  to  be  on  the  higher  side, 
recording  a  velocity  higher  than  that  expected  to  exist  at 
that  distance  from  the  wall.  Such  a  shift  has  been  observed 
also  by  other  workers and  is  suggested  as  being  due 
to  a  high  velocity  gradient.  A  mathematical  illustration  has 
been  given  by  Knudsen  and  Katz^^^^ .  A  correction  of  0.15  -  0.18 
of  the  leading  dimension  of  the  probe  tip  has  been  recommended 
by  Livesey^^^^.  On  applying  this  correction  to  the  measured 
distance,  i.e.  adding  0.15  of  the  leading  dimension  of  the 
probe  to  the  measured  distance  and  assuming  the  observed 
velocity  to  exist  at  this  corrected  distance  from  the  wall, 
the  observations  were  again  in  satisfactory  agreement  with 
Nikuradse 's  correlation. 

For  all  the  calculations  the  plastic  tube  was  considered 
to  be  a  smooth  tube  and  for  the  determination  of  the  friction 
velocities  the  usual  Moody  chart  was  used. 

No  observations  were  made  for  laminar  flow,  Reynolds 
number  less  than  2300,  because  the  length  of  the  tube  used 
was  insufficient  for  fully  developed  laminar  flow  and  also 
because,  since  the  main  area  of  interest  was  the  turbulent 
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region,  such  a  setup  for  laminar  flov7  was  considered 
unnecessary. 

Computer  Programme; 

The  theoretical  approach  of  Deissler  and  Taylor 
was  applied  to  the  rectangular  channel.  Instead  of  using  the 
graphical  iteration  proposed  by  the  authors,  an  attempt  was 
made  to  write  a  computer  programme  which  would  iterate  in  the 
same  way.  As  a  first  approximation  the  universal  velocity 
distribution  correlation  proposed  by  Deissler  was  represented 
by  simple  semilogarithmic  equations.  In  the  report  of  Deissler 
and  Taylor,  certain  new  dimensionless  parameters  have  been 
suggested  for  simplifying  the  computation  work.  However, 
since  the  work  was  to  be  done  by  a  computer,  these  new  para- 
mete3:s  offered  no  advantage  and  hence  were  not  used.  The 
method  involves  the  computation  of  the  wall  shear  stress  by 
a  force  balance  between  successive  velocity  gradient  lines. 

The  areas  can  be  computed  quite  easily  if  equations  are  fitted 
to  the  velocity  gradient  lines.  The  prograirnne  v/as  initially 
written  in  FORGO  for  the  IBM  1620  digital  computer  at  the 
University  of  Alberta.  However,  the  memory  of  the  computer 
was  insufficient  for  the  computer  to  be  able  to  read-in  the 
entire  programme  when  written  in  FORGO.  The  programme  v/as 
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then  re-written  in  FORTRAN  II  for  the  IBM  1620.  A  description 
of  the  programme  is  given  in  Appendix  II. 

A  computer  running  time  of  approximately  45  minutes 
was  required  to  go  through  3  iterations.  The  programme  was 
not  run  to  its  final  result  because  of  a  lack  of  available 
computer  time.  Assuming  that  about  20  iterations  may  be 
required  to  reach  the  values  within  the  margin  of  error  allowed, 
a  running  time  of  about  5  to  6  hours  will  be  necessary  on 
the  IBM  1620  computer. 

A  major  difference  between  the  method  of  Deissler  and 
Taylor  and  that  used  to  write  the  computer  programme  may  be 
mentioned  here.  Deissler  and  Taylor  assume  that  "expressions 
for  eddy  diffusivity  that  had  been  verified  for  flow  and  heat 
transfer  in  tubes  apply  in  general  along  lines  normal  to  a 
wall".  However,  it  was  suggested  that  the  expressions  need 
not  apply  along  lines  normal  to  a  wall  in  non- circular 
passages  and  that  a  better  assumption  would  be  the  applicability 
of  the  expressions  along  a  velocity  gradient  line.  This  latter 
suggestion  was  used  in  the  computer  programme.  A  more  complete 
run  of  the  programme  is  necessary  before  the  superiority  of 
this  assumption  can  be  verified. 

Another  major  point  of  objection  to  the  method  of 
Deissler  and  Taylor  is  its  inability  to  recognize  and  to 
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account  for  any  secondary  flow.  Deissler  and  Taylor  have 
accepted  this  limitation  of  their  method  but  have  also  pointed 
out  the  lack  of  evidence  of  secondary  flow  in  the  form  of 
experimental  data  except  that  presented  by  Nikuradse  . 

The  experimental  results  on  the  channel  in  the  present  investi¬ 
gation  gave  curves  similar  to  those  of  Nikuradse  with  the 
constant  velocity  lines  showing  a  tendency  to  be  pushed  towards 
the  corners  and  away  from  the  walls.  According  to  the  method 
of  Deissler  and  Taylor,  the  velocity  distribution  at  the 
centre  line  along  the  width  of  the  channel  cross-section 
would  be  that  given  by  the  universal  velocity  distribution 
for  parallel  plates.  This  theoretical  curve  was  calculated 
for  the  same  maximum  velocity  in  the  traverse,  as  described 
in  the  next  section.  The  maximum  deviation  of  this  curve 
from  the  experimentally  determined  curve  would  be  a  good 
indication  of  the  effect  of  the  secondary  flow.  This  deviation 
for  a  Reynolds  number  of  13,650  was  found  to  be  +17.4%  of  the 
theoretical  value. 

Because  of  the  incomplete  success  of  the  computer 
programme,  the  graphical  approach  was  also  used.  Computations 
were  made  for  one  arbitrary  value  of  the  pressure  gradient  and 
the  result  is  shown  in  Figure  4.  The  procedure  appears  to  be 
rather  insensitive  to  the  shape  of  the  assumed  velocity 
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gradient  lines  since  the  value  of  the  calculated  shear  stress 
depends  only  on  the  ratio  of  the  area  between  two  successive 
velocity  gradient  lines  and  the  distance  between  them  at  the 
wall. 

The  results  of  the  graphical  computation  do  not  show 
any  tendency  for  the  constant  velocity  lines  to  shift  towards 
the  corners.  Thus  it  is  felt  that  the  method  proposed  by 
Deissler  and  Taylor  probably  does  not  take  care  of  all  the 
factors  influencing  the  flow  in  a  passage  of  non-circular 
cross-section. 

The  computer  programme  was  temporarily  discontinued 
until  more  information  about  secondary  flow  could  be  obtained. 

It  may  be  possible  to  compare  the  experimental  results  which 
indicate  the  presence  of  secondary  flow  with  the  constant 
velocity  curves  obtained  from  the  computer  programme.  Such 
a  comparison  may  suggest  a  method  of  characterization  of  the 
secondary  flow. 

In  the  method  of  Deissler  and  Taylor,  the  authors  require 
lines  of  symmetry  to  be  drawn  initially.  Any  such  line  of 
symmetry  is  accepted  to  be  a  velocity  gradient  line  and  the 
trial  and  error  procedure  of  iteration  to  get  better  velocity 
gradient  lines  is  applied  to  velocity  gradient  lines  drawn  by 
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the  eye  in  between  two  lines  of  symmetry.  In  the  case  of  the 
computer  programme,  however,  the  diagonals  of  the  cross  section 
were  assumed  to  be  velocity  gradient  lines.  In  the  light  of 
Prandtl's  suggestion  that  the  secondary  flow  in  a  channel 
occurs  along  the  lines  bisecting  the  corners,  the  velocity 
gradient  lines  may  be  assumed  to  be  along  the  bisectors 
instead  of  the  diagonals.  Such  an  assumption  may  give  a  better 
picture  of  the  velocity  distribution  in  an  actual  channel, 
although  it  does  not  take  into  account  the  effect  of  the 
secondary  flow. 

Channel  Experiments ; 

Velocity  distributions  were  determined  at  Reynolds 
numbers  of  5070,  11,700,  13,650  and  20,500,  The  data  from  the 
vertical  traverses  made  at  a  Reynolds  number  of  13,650  have 
been  plotted  in  Figures  5  to  11,  Certain  constant  velocity 
lines  across  the  cross  section  have  also  been  plotted  from 
these  figures.  Velocity  distributions  at  the  other  Reynolds 
numbers  follow  a  similar  pattern. 

An  examination  of  the  velocity  distribution  at  the 
same  point  of  traverse  for  different  Reynolds  numbers  shows 
that  the  effect  of  the  secondary  flow  probably  moves  towards 
the  centre  of  the  channel  with  increasing  Reynolds  numbers. 
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This  is  indicated  by  the  shifting  of  the  points  of  inflexion 
on  the  velocity  distribution  curves. 

The  shape  of  these  velocity  distribution  curves,  with 
their  points  of  inflexion,  is  not  incompatible  with  the  nature 
of  the  secondary  flow  proposed  by  Prandtl  after  a  visual  study. 
A  qualitative  cross-sectional  diagram  of  the  flow  indicating 
secondary  flow  can  be  used  to  project  a  velocity  distribution 
curve.  Figure  12  is  an  illustration  shov/ing  the  agreement  of 
Prandtl 's  picture  of  secondary  flow  and  the  velocity  distri¬ 
bution  along  a  vertical  line  on  the  cross  section  of  the 
conduit.  The  curve  demonstrates  a  qualitative  agreement 
with  the  secondary  flow  indicated  by  the  arrows.  In  the  region 
close  to  the  corners  the  secondary  flow  causes  fluid  particles 
of  a  higher  momentum  to  travel  along  the  bisector  of  the 
corner  towards  the  corner;  while  along  the  center  line  fluid 
particles  of  lower  momentum  rise  up  to  the  centre.  Thus  on 
a  certain  vertical  traverse  one  may  expect  particles  at  certain 
points  to  have  velocities  higher  than  those  which  would  exist 
in  the  absence  of  such  a  transfer  of  momentum;  and  at  points 
the  velocities  would  be  lower  than  those  expected  in  the 


absence  of  secondary  flow. 
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It  is  obvious  from  the  velocity  distribution  curves 
that,  in  spite  of  the  favoured  transfer  of  momentum  towards 
the  corner,  at  no  stage  does  the  velocity  nearer  the  wall 
exceed  that  at  a  greater  distance;  so  that  there  is  only  one 
maximum  for  the  curve,  namely  at  the  centre. 

In  Figures  13  and  14,  some  of  the  experimental  data 
have  been  plotted  on  a  dimensionless  u”*"  versus  y*"  co-ordinate 
system.  To  convert  the  data  to  these  co-ordinates  it  is 
necessary  to  calculate  the  friction  velocity  u*.  No'w, 

U*  = 

so  that  a  knowledge  of  the  shear  stress  distribution  along 
the  wall  is  essential.  Unlike  as  in  a  circular  conduit,  the 
shear  stress  along  the  wall  in  non-circular  conduits  is  not 
the  same  everywhere.  In  an  attempt  to  determine  the  extent 
of  variation  of  the  shear  stress  along  the  wall,  in  Figure  15 
u/ujYiax  been  plotted  against  y/b  for  two  traverses  where 
u  is  the  point  velocity,  Uj^^x  maximum  velocity  of  the 

traverse,  y  is  the  distance  from  the  wall  and  b  is  half  the 
height  of  the  channel.  Although  no  observations  for  the 
laminar  sublayer  have  been  made,  the  observed  data  may  be 
extrapolated  to  obtain  a  rough  indication  of  the  variation  of 
the  shear  stress.  In  the  laminar  region  the  wall  shear  stress. 
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Tq,  is  connected  to  the  velocity  gradient  by  Newton's  law 


of  viscosity 


(tq)  1  ^  (du/dy)  1 

(Tq)  2  (du/dy)  2 

Since  in  Figure  15,  has  been  plotted  against  y/b 


slope 

-  du 

X  ^ 

dy 

^ax 

du  _ 

slope 

,,  ^max 

dy 

b 

Hence 

(”^0^  1  „  (slope)  ^  ^  (^max)  1 

(xq)  2  (slope)  2  (u^ax)  2 

Comparing  the  centre-line  shear  stress  with  that  at  3.75 
inches  from  the  centre,  from  Figure  15 


("0)2  =  <"0)1 

Thus  the  wall  shear  stress  3.75  inches  from  the  centre  is  87% 
of  that  at  the  centre.  However,  since  the  calculations  are 
based  on  extrapolated  data,  this  may  be  taken  only  as  an 
indication  of  the  trend  of  the  behaviour  of  the  wall  shear 
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Apparently,  it  is  not  possible  to  calculate  the  actual 
wall  shear  stress  at  a  point  from  the  overall  pressure  drop. 

An  average  value,  if  used,  is  likely  to  give  a  distorted  picture 
of  the  velocity  distribution  on  the  u"*"  versus  y"*"  co-ordinates. 
However,  one  may  assume  Deissler's  universal  velocity  distri¬ 
bution  law  for  the  turbulent  region  to  apply  and,  combining 
this  with  the  known  value  of  the  velocity  at  the  centre  of  a 
certain  traverse,  it  is  possible  to  calculate  by  trial  and 
error  the  friction  velocity  at  the  wall  for  that  traverse. 

This  procedure  of  calculation  has  been  used  for  the  plots  in 
Figures  13  and  14.  Since  the  effect  of  the  secondary  flow 
will  be  the  least  on  the  traverse  at  the  centre,  the  data 
may  be  expected  to  agree  with  the  universal  velocity  distri¬ 
bution.  Such  is  the  case  as  can  be  seen  from  Figure  13. 

For  any  other  traverse  the  unequal  distances  from  the  side 
walls  makes  the  method,  of  calculating  u*  using  the  Deissler 
correlation  for  infinite  parallel  plates,  questionable. 

Figure  14  shows  the  deviation  of  the  observed  data  from  the 
correlation  for  parallel  plates,  as  expected. 

However,  a  plot  such  as  Figure  13  may  be  used  for 
studying  the  effect  of  Reynolds  number  on  the  velocity  distri¬ 
bution.  The  "average  slope"  of  the  curves  decreases  with 
decreasing  Reynolds  numbers. 
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In  Figure  13,  although  the  observed  points  seem  to  be 
in  reasonable  agreement  with  Deissler's  universal  velocity 
distribution,  the  non-linearity  of  the  data  is  apparent. 

This  suggests  that  the  effect  of  the  side  walls  must  be  con¬ 
sidered  in  channel  flow.  in  the  present  set  of  experi¬ 
ments  all  the  data  were  taken  for  an  aspect  ratio  of  8:1. 
Determination  of  the  influence  of  the  aspect  ratio  on  the 
secondary  flow  induced  in  a  channel  would  require  a  study  with 
a  range  of  aspect  ratios  varying  from  1:1  to  the  case  where 
an  assumption  of  infinite  parallel  plates  is  valid,  Nikuradse's 
data  are  for  a  rectangular  channel  0.8  cm.  x  2.8  aa.  in  cross 
section  i.e.  an  aspect  ratio  of  3.5:1.  However,  no  method 
of  characterizing  the  secondary  flow  has  been  suggested.  Such 
a  characteristic  quantity  pertaining  to  the  secondary  flow 
is  likely  to  involve  the  aspect  ratio  and  also  the  Reynolds 
number;  both  of  which  seem  to  influence  the  magnitude  and 
penetration  of  the  secondary  flow. 

An  examination  of  the  friction  factor  data  of  Charles 
on  the  same  channel,  reveals  a  slight  deviation  from  the 
Blasius  equation 

f  =  0.0792  (Re)"^"^^' 

The  difference  is  greater  at  lower  Reynolds  numbers.  An 
attempt  can  be  made  to  explain  this  on  the  basis  of  the  secon¬ 
dary  flow.  As  observed  earlier,  the  average  slope  of  the 
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data  on  the  versus  plot  is  lower  for  lower  Reynolds 
numbers,  with  the  data  approaching  Deissler's  curve  at  higher 
Reynolds  numbers.  It  may  be  suggested  that  the  magnitude  of 
the  secondary  flow,  is  comparable  to  the  magnitude  of  the 
flow  in  the  main  direction  at  lower  Reynolds  numbers.  As 
the  Reynolds  number  increases  the  relative  magnitude  and  hence, 
the  effect  on  the  main  flow  diminishes.  This,  consequently, 
results  in  a  pressure  drop  nearly  the  same  as  in  the  flow 
of  a  fluid  through  a  tube  of  equivalent  diameter,  there  being 
no  secondary  flow  in  tubes.  Thus  the  pressure  drop  due  to 
the  secondary  flow  seem.s  to  be  significant  when  compared  with 
the  pressure  drop  caused  by  the  main  flow,  at  low  Reynolds 
numbers  in  the  turbulent  region.  At  higher  Reynolds  numbers, 
however,  the  pressure  drop  due  to  the  secondary  flow  forms  a 
less  significant  part  of  the  total  pressure  drop.  The  result 
is  a  higher  friction  factor  than  that  for  tubes  at  the  lower 
Reynolds  numbers,  the  difference  gradually  decreasing  with 
increasing  Reynolds  numbers. 
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CONCLUSIONS  AND  RECOMMENDATIONS 

1.  An  attempt  to  correlate  the  observed  velocity  distribution 
for  the  rectangular  conduit  with  u"^  -  y’*'  parameters 
raises  the  difficulty  of  calculating  the  correct  friction 
velocity,  u*.  Such  an  attempt  in  this  work  was  unsuccess¬ 
ful. 

2.  Secondary  flow  was  observed  to  exist  in  turbulent  flow 
through  the  rectangular  conduit  of  an  aspect  ratio  8:1. 

A  shift  of  the  constant  velocity  lines  towards  the  corners 
was  observed,  which  suggests  such  an  existence  of  secondary 
flow. 

3.  The  velocity  distribution  does  not  seem  to  obey  the 
universal  velocity  distribution  law.  The  vertical  tra¬ 
verses  revealed  points  of  inflexion  in  the  velocity  dis¬ 
tribution  curves  which  cannot  be  accounted  for  by  a 
straight  line  correlation  on  a  semi- logarithmic  plot. 

4.  The  inflexions  in  the  velocity  distribution  curves  of  the 
vertical  traverses  are  not  incompatible  with  Prandtl's 
suggestion  that  the  secondary  flow  causes  fluid  particles 
of  higher  momentum  to  travel  towards  the  corners  along 
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the  diagonals. 

5.  The  u"*"  versus  y**"  plot  of  the  centre  line  data  suggests 
that  the  magnitude  of  the  secondary  flow  relative  to  the 
mean  flow  decreases  with  increasing  Reynolds  number,  since 
the  deviation  of  the  curves  at  lower  Reynolds  numbers  is 
greater  from  the  curve  for  parallel  plates. 

6.  The  Deissler  and  Taylor  approach  for  the  theoretical  cal¬ 
culation  of  the  turbulent  velocity  distribution  in  non¬ 
circular  conduits  does  not  consider  secondary  flow. 

Hence  it  probably  does  not  represent  the  actual  flow 
mechanism  in  conduits  of  non-circular  cross  section, 

7.  A  general  velocity  distribution  correlation  for  the 
channel  may  be  required  before  a  characterization  para¬ 
meter  for  the  secondary  flow  can  be  detemined, 

8.  A  characteristic  quantity  for  the  secondary  flow  is 
expected  to  involve  the  aspect  ratio,  and  also,  the 
Reynolds  number. 

9.  Conclusion  1,  immediately  suggests  that  investigation  be 
carried  out  for  direct  measurements  of  the  wall  shear 
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10.  From  conclusion  3,  it  may  be  suggested  that  more  data 

be  obtained  at  Reynolds  numbers  other  than  those  studied 
in  this  vvork. 

11.  Further  study  with  various  aspect  ratios  is  recommended. 
It  will  be  of  interest  to  determine  the  minimum  aspect 
ratio  at  which  the  turbulent  velocity  distribution  at  the 
centre  line  of  the  cross  section  is  the  same  as  that  for 
infinite  parallel  plates. 
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APPENDICES 

I.  Calibration  of  the  Pitot  Tube; 

The  preliminary  experiments  of  the  project  were  designed 
to  determine  the  adequacy  of  the  pitot  tube  as  a  mean  velocity 
measuring  instrument.  Flow  of  fluids  through  conduits  of 
circular  cross  section  has  been  studied  extensively  and  the 
mean  velocity  distribution  in  them  is  well  understood.  Experi¬ 
ments  were,  therefore,  conducted  with  a  tube  of  dimensions 
of  the  same  order  as  those  of  the  channel. 

The  deteinnination  of  the  velocity  distribution  and 
subsequent  integration  over  the  cross  section  can  be  used  to 
check  the  adequacy  and  determine  the  coefficient  of  the  pitot 
tube,  respectively.  For  a  pitot  tube  with  a  smooth  conical 
head  and  a  diameter  ratio  of  the  magnitude  as  in  the  present 
case,  the  coefficient  has  been  reported  to  be  unity . 

The  experimental  set-up  consisted  of  a  constant  level 
tank  on  a  ten- foot  structure  discharging  at  a  constant  flow 
rate  into  a  vertical,  6  foot  long,  0,744  inch  inside  diameter, 
plastic  tube.  The  pitot  tube  of  0.018  inches  outside  diaiifieter, 
was  introduced  4-1/2  feet  from  the  upstream  end  of  the  tube, 

A  length  of  72  diameters  was  considered  satisfactory  enough 
for  achieving  fully  developed  flow.  Such  a  fully  developed 
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flow  is  necessary  to  be  able  to  determine  how  well  the  pitot 
tube  registers  the  velocity  at  different  distances  from  the 
wall;  which  would  be  a  test  for  the  adequacy  of  the  pitot 
tube  as  a  mean  velocity  measuring  instrument.  For  straightening 
the  flow  and  inducing  turbulence,  a  bundle  of  glass  tubes  was 
fixed  at  the  entrance  to  the  plastic  tube.  Flow  rate  was 
controlled  by  a  globe  valve  and  measured  by  a  rotameter.  The 
rotameter  calibration  of  Wood was  used  after  checking  it 
with  a  tank  on  a  weigh-scale.  The  pressure  difference  read¬ 
ings  v/ere  taken  with  a  manometer- cathetome ter  set-up. 

Experimental  data  on  the  tube  have  been  tabulated  in 
Appendix  III. 

II.  Computer  Programme : 

The  computer  programme  designed  to  calculate  the 
theoretical  velocity  profiles  based  on  the  method  suggested 
by  Deissler  and  Taylor has  been  described  in  this  section. 

It  was  initially  attempted  in  FORGO,  but  due  to  the  lack  of 
sufficient  memory  in  the  IBM  1620  computer  at  the  University 
of  Alberta  it  was  re-written  in  FORTRAN  II  with  which  the 
computer  has  a  longer  memory  available  for  use. 

The  step-wise  procedure  of  Deissler  and  Taylor  is  as 


follows : 
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(i)  The  lines  of  synimetry  are  drawn 

(ii)  Assumed  velocity  gradient  lines  are  drawn 

(iii)  From  the  velocity  gradient  lines  the  constant 
velocity  lines  are  calculated. 

(iv)  A  new  and  more  accurate  set  of  velocity  gradient 
lines  is  then  drawn  normal  to  the  constant 
velocity  lines.  The  method  thus  requires  a 
trial  and  error  procedure, 

A  force  balance  for  an  element,  as  shown  in  Figure 
gives  the  following  relation  for  the  wall  shear  stress, 

_  „  /  dP  X 

AL  ^  dz  ^ 
dP 

where  is  the  pressure  gradient  uniform  over  the  cross 
dz 

section  and  the  mean  flov/  is  along  the  z  axis. 

The  computer  programme  goes  through  the  following 

steps : 

(i)  The  initial  data  comprising  the  points  on  the 

assumed  velocity  gradient  lines  are  read  in.  Values  of 

~  are  also  read  in. 
dz 

(ii)  The  choice  of  the  form  of  equations  fitted  to  the 
initial  and  subsequently  generated  data,  was  governed  by 
simplicity  and  adequacy  considerations.  A  linear  form  obviously 
is  inadequate.  A  second  order  equation  may  not  be  adequate 
in  that  it  may  cause  over-smoothening  of  the  data.  A  third 
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order  fit  was,  therefore,  tried.  The  equations  fitted 
were  of  the  form 

1  2 

y  =  Ax  +Bx  +Cx+D 

and  the  least  squares  criterion  was  used  in  the  procedure. 

The  x-axis  was  taken  as  co- incident  with  the  lower  side  of 
the  rectangular  cross-section  and  the  y-axis  was  the  vertical 
at  the  centre  of  the  cross-section. 

(iii)  Areas  between  successive  velocity  gradient  lines 
are  calculated. 

(iv)  The  values  of  (iii)  are  used  for  calculating  Tq, 
u'‘'max'  values  as  certain  fractions  of  the  calculated 
and  the  y"*"  and  the  y  values  for  the  calculated  u"^  values. 

For  the  calculation  of  u'’‘max'  Deissler  and  Taylor  univer¬ 

sal  velocity  distribution  correlation  is  used. 

(v)  The  X  values  for  the  y  values  determined  in  (iv) 
require  the  solution  of  the  equations  of  the  two  velocity 
gradient  lines,  the  area  betv/een  which  has  been  used  for 
calculating  the  y  value.  The  solution  is  achieved  by  Newton's 
method  of  successive  approximation.  The  x  value  used  is  the 
arithmetic  mean  of  the  two  solutions, 

(vi)  Equations  are  fitted  to  the  x,  y  values  in  the 


same  manner  as  in  (ii) , 
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(vii)  Points  are  arbitrarily  chosen  on  the  constant 
velocity  line  closest  to  the  centre  and  points  on  successively 
lower  lines  are  determined  at  which  the  product  of  the  slope 
of  the  line  from  the  point  on  the  preceding  constant  velocity 
line  and  the  slope  of  the  constant  velocity  line  is  -1 
+0.05.  The  tolerance  of  j:0.05  has  been  chosen  quite  arbi¬ 
trarily. 

(viii)  Equations  are  fitted  to  these  points  as  in  (ii) , 


These  steps  are  repeated  until  in  tv/o  successive  compu¬ 


tations  the  co-ordinates  of  the  points  of  constant  velocity 
do  not  differ  by  more  than  a  preset  figure.  The  choice  of 
this  figure  depends  on  the  accuracy  sought  and  the  size  of 
the  channel.  In  the  present  instance  a  figure  of  0.005 
inches  was  chosen  for  the  1  inch  high  channel, 

III  Data  with  Tube  -  Preliminary  Experiments; 


Length  of  the  tube 


6  ft. 


Inside  diameter  of  the  tube 


0,744  in. 


Outside  diameter  of  the  pitot  tube 


0.018  in. 


Inside  diameter  of  the  pitot  tube 


0,009  in. 


Entrance  length  of  the  tube 


4-1/2  ft. 


Temperature  of  the  water 


68°F 


Assumed  pitot  tube  coefficient 


1,000 
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Table  1 


Flow  rate  of  water  =  0.1007  Ibs/sec 


Reynolds  number  =  2700 


Distance  from 

Pressure  Diff. 

Velocity 

u+ 

Tube  Wall 

h  cm  of  water 

u  ft/sec 

yu* 

u 

^in. 

u* 

0.0342 

0.110 

0.482 

9.2 

12.17 

0.2534 

0.185 

0.606 

68.25 

15.8 

0.3720 

0.190 

0.623 

100 

16 

Distance 
Tube  Wall 

yin. 

Table 

Flow  rate  of  water 

Reynolds  number 

from  Pressure  Diff. 

h  cm  of  water 

2 

=  0,2692  Ibs/sec 

=  9150 

Velocity  y"^ 

u  ft/sec  yu* 

u-' 

u 

u* 

0.0257 

0.880 

1.278 

19,91 

15,25 

0.0408 

0.835 

1.235 

31.62 

14.86 

0.0500 

0.855 

1.250 

38,75 

15.03 

0.0500 

0.910 

1.29 

38.75 

15.5 

0.0750 

0,980 

1.339 

58,1 

16.1 

0.0750 

0.980 

1.339 

58,1 

16.1 

0.1250 

1.130 

1.438 

96.9 

17.3 

0.1250 

1.120 

1.433 

96,9 

17.13 

0.1710 

1.250 

1.512 

135.8 

13. 18 

0.1750 

1,250 

1.512 

135,8 

18,18 

0.2249 

1.335 

1,592 

174.5 

19,15 

0.2250 

1.370 

1.583 

174,5 

19,05 

0.2750 

1 . 460 

1.635 

213,0 

19.65 

0.2750 

1.440 

1.623 

213 

19.5 

0.3250 

1.515 

1.663 

252 

20.0 

0.3250 

1.510 

1.662 

252 

20,0 

0.3720 

1.540 

1.683 

291 

20.21 

-  ^-.4  - 


Table  3 


Flow  rate  of  water  =  0.52  Ibs/sec 
Reynolds  number  =  17690 


Distance  from 
Tube  Wall 

Yin. 

Pressure  Diff. 
h  cm  of  water 

Velocity 
u  ft/sec 

-4 

u-*- 

u 

u* 

0.0296 

3.455 

2.701 

45.05 

16.95 

0.0314 

3.175 

2.590 

47.5 

16.26 

0.0456 

3.290 

2,638 

67,1 

16.55 

0.0606 

3 . 480 

2,714 

38.0 

17.02 

0.0682 

3.915 

2,878 

98.4 

18.05 

0.0790 

3.705 

2.800 

113.3 

17.18 

0.0949 

4.200 

2.930 

135.3 

18,7 

0.1040 

4.025 

2.92 

148.5 

18.31 

0.1230 

4.250 

2.998 

181.9 

18.8 

0.1347 

4.580 

3.122 

190.2 

19,53 

0.1580 

4.565 

3.106 

223 

19.49 

0.1352 

5,005 

3.254 

260,5 

20.4 

0.1930 

4.850 

3.202 

271 

20,1 

0.2245 

5.280 

3,342 

314.4 

20.98 

0.2200 

5,100 

3.282 

308 

20.6 

0.2480 

5.280 

3,343 

347 

20.93 

0.2650 

5,525 

3.418 

370.5 

21.44 

0.2746 

5.440 

3.398 

383.6 

21.3 

0.2967 

5.565 

3,426 

414 

21.5 

0.2943 

5.630 

3.450 

411 

21.64 

0.3098 

5 , 640 

3.455 

432.5 

21.68 

0.3199 

5.720 

3.482 

446,5 

21.83 

0.3272 

5.690 

3.472 

456.5 

21.78 

0.3449 

5.770 

3.490 

477 

21.88 

0.3507 

5.760 

3.486 

489 

21.86 

0.3703 

5.770 

3.490 

516.5 

21.88 

0.3720 

5.775 

3,493 

518 

21.90 
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Table  4 


Flow  rate  of  water  =  0,65  lb. /sec. 
Reynolds  number  =  22260 


Distance  from 
Tube  Wall 

Yin. 

Pressure  Diff. 
h  cm  of  water 

Velocity 
u  ft/sec 

yu* 

Sv 

u-^ 

u 

u* 

0.0302 

4.475 

3.078 

51.3 

15.69 

0.1046 

5.995 

3.564 

117.8 

18.17 

0.1790 

6.975 

3.843 

304.3 

19.6 

0.2534 

7.810 

4.070 

430.5 

20.73 

0.3278 

8.200 

4.170 

557,0 

21,22 

0.3720 

8,380 

4.210 

632,0 

21,46 

Table  5 


Flow  rate  of  water  =  1.00  lb. /sec, 
Reynolds  number  =  34,300 


Distance  from 
Tube  Wall 

Yin. 

Pressure  Diff. 
h  cm  of  water 

Velocity 
u  ft/sec 

yu* 

-0 

u+ 

u 

u* 

0.009 

10.075 

3.576 

22.4 

12.5 

0.1488 

25.770 

5.725 

370 

20,0 

0.2232 

29.055 

6,075 

555 

21.25 

0.2976 

32,335 

0.410 

740 

22.40 

Table  6 


Average  Velocity 
(actual) 
in/sec 


Average  Velocity  %  Error 

(integrated) 

in/sec  


17,14 

33.05 

41.4 


17.31  +  0.992 
33.05  0.000 
40.93  -  1,139 
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IV  Data  with  the  Channel; 

Inner  dimensions  of  the  channel  at 

the  test  cross-section:  Width  =  8.000  inches 

Height  =  1.000  inches 

Outside  diameter  of  the  pitot  tube  =  0.016  inches 

Inside  diameter  of  the  pitot  tube  =  0.008  inches 

Pitot  tube  coefficient  =  1.000 

Inside  diameter  of  the  static  top  =  0.0625  inches 

Specific  gravity  of  manometric  fluid  =  1.0895 

Table  7 

Flow  rate  of  water  =  7.3  USGPM 

Mean  temperature  of  water  =  75 ®F  (j;2®F) 

Reynolds  number  (based  on  equivalent  diameter)  =  5070 

Traverse  No.  1; 

Distance  of  the  vertical  traverse  from  the  centre  =  0.000  in. 


Distance  from 
Upper  Wall  ,  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0080 

0,641 

0.3482 

0.0375 

0.839 

0.3973 

0.0750 

0.989 

0.4320 

0.1750 

1.130 

0.4620 

0.2750 

1.170 

0.4700 

0.3750 

1.262 

0.4885 

0.5000 

1.286 

0.4930 

0.6250 

1.262 

0.4885 

0.7250 

1.170 

0.4700 

0.8250 

1.130 

0.4620 

0.9250 

0.989 

0.4320 
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Traverse  No.  2; 


Distance  of  the  vertical  traverse  from  the  centre  =  1.000  in. 


Distance  from 
Upper  Wall,  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0080 

0.0250 

0.700 

0.3313 

0.0500 

0.836 

0.3976 

0.0750 

0.860 

0.4030 

0.1750 

0.975 

0,4290 

0.2750 

1.025 

0.4405 

0.3950 

1.066 

0 . 449 

0.5000 

1.089 

0.453 

0.6250 

1.066 

0 . 449 

0.7250 

1.024 

0.4405 

0.8250 

0.975 

0,4290 

0.9250 

0.860 

0,4030 

Traverse  No.  3: 

Distance  of  the  vertical  traverse  from  the  centre  =  2.000  in. 


0.0080 

0,0250 

0,411 

0.2790 

0.0500 

0.493 

0,3050 

0.0750 

0,615 

0,3410 

0.1750 

0.687 

0.3592 

0.2750 

0.810 

0.3913 

0.3750 

0.945 

0.4225 

0.5000 

0.995 

0,4340 

0.6250 

0.952 

0.4240 

0.7250 

0.810 

0,3913 

0.8250 

0.687 

0.3592 

0.9250 

0.613 

0.3410 
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Traverse  No.  4 

Distance  of  the  vertical  traverse  from  the  centre  =  3.000  in 


Distance  from 
upper  Wall,  in. 


Pressure  Diff. 
Ah  cm  of  fluid 


Velocity 
u  ft/sec 


0.0080 

0.0250 


Traverse  No.  5 


0.0500 

0.198 

0.1935 

0.0750 

0.346 

0.2556 

0.1750 

0.455 

0.2935 

0.2750 

0.635 

0.3460 

0.3750 

0.672 

0.3570 

0.5000 

0,712 

0.3638 

0.6250 

0.672 

0.3570 

0.7250 

0.635 

0.3460 

0.8250 

0.455 

0.2935 

0.9250 

0.346 

0.2556 

the  vertical 

traverse  from  the 

centre  =  3. 

0.0080 

0.0250 

0.060 

0.1058 

0.0500 

0.224 

0.2060 

0,0750 

0.326 

0.2484 

0-1750 

0.477 

0.3000 

0.2750 

0.553 

0,3235 

0.3750 

0.598 

0.3362 

0.5000 

0.630 

0.3452 

0.6250 

0.598 

0,3362 

0.7250 

0.553 

0.3235 

0.8250 

0.477 

0.3000 

0.9250 

0.326 

0.2484 
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Traverse  No.  6 


Distance  of  the  vertical  traverse  from  the  centre  =  3.750  in. 


Distance  from 
Upper  Walh  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0030 

0.0250 

0.0500 

0.050 

0.0971 

0.0750 

0.126 

0.1545 

0.1750 

0.211 

0.2000 

0.2750 

0.236 

0.2110 

0.3750 

0,253 

0.219 

0.5000 

0.253 

0.219 

0.6250 

0.253 

0,219 

0.7250 

0,236 

0.2110 

0.8250 

0.211 

0.2000 

0.9250 

0.126 

0.1545 

„  0002. 0 
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Table  8 

Flow  rate  of  water  =  18  USGPM 

Mean  temperature  of  water  =  80 ®F  (±2°F) 

Reynolds  number  (based  on  equivalent  diameter)  =  11,700 


Traverse  No.  1 


Distance  of  the  vertical  traverse  from  the  centre  =  0.809  in. 


Distance  from 
Upper  Wall 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0,0253 

1,820 

0,582 

0,0535 

3,065 

0,756 

0.0785 

3.570 

0,816 

0.1035 

4,010 

0.864 

0.1535 

4.980 

0.964 

0.2035 

5.590 

1.021 

0,3035 

6.365 

1.090 

0.4035 

6.700 

1.118 

0.5000 

7,045 

1,144 

0,6035 

6,665 

1,115 

0.3035 

5.590 

1,021 

0.9035 

4,075 

0,872 

0.9855 

2.055 

0,619 
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Traverse  No.  2 


Distance  of  the  vertical  traverse  from  the  centre  =  0.9091  in 

Distance  from  Pressure  Diff.  Velocity 

Upper  Wall,  in.  Ah  cm  of  fluid  u  ft/sec 


0.0199 

2.020 

0.614 

0.0481 

2.965 

0.744 

0.0731 

3.670 

0.828 

0.0981 

4.040 

0.869 

0.1981 

3.245 

0.950 

0,2981 

3.555 

1.000 

0.3981 

3.870 

1,041 

0.4981 

5.895 

1,054 

0.5000 

5.915 

1,058 

0.5981 

5.730 

1,040 

0.6981 

5.250 

0,995 

0,7981 

4.835 

0.956 

0.8981 

4.310 

0.896 

0.9801 

2.190 

0,639 

Traverse  No»  3 

Distance  of  the  vertical  traverse  from  the  centre  =  1,9091  in 


0,0217 

1,530 

0,538 

0.0481 

2.090 

0.628 

0.0731 

2.690 

0.712 

0.0981 

3,210 

0,779 

0.1481 

3.755 

0 . 842 

0.1981 

4.090 

0.879 

0.2981 

4.990 

0.971 

0.3981 

5,380 

1.009 

0.4981 

5.365 

1.007 

0.5000 

5.375 

1.008 

0.5981 

5.275 

1.000 

0.5981 

4.870 

0.960 

0,7981 

4,410^ 

0.885 

0.8981 

3.395 

0.8005 

0.9801 

1.375 

0.510 
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Traverse  No.  4 


Distance  of  the  vertical  traverse  from  the  centre  =  2.9091  in. 


Distance  from 
Upper  Wall,  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0191 

1.275 

0.4905 

0.0443 

1.285 

0.493 

0.0693 

1.705 

0.567 

0.0943 

2.200 

0.645 

0.1443 

3.190 

0.776 

0.1943 

3.535 

0.816 

0.2943 

4.385 

0.910 

0.3943 

4.775 

0.950 

0.5000 

5.305 

1.001 

0.5943 

4.815 

0,954 

0.6943 

4.385 

0.910 

0.7943 

3.720 

0,838 

0.8943 

2.355 

0,667 

0.9821 

1.030 

0.441 

Traverse  No.  5 

Distance  of  the  vertical 

traverse  from  the 

centre  =  3.‘ 

0.0693 

2.110 

0,631 

0.0943 

2.460 

0.681 

0.1943 

3.070 

0.762 

0.2943 

3.810 

0.849 

0.3943 

4.410 

0.913 

0.5000 

4.600 

0.9325 

0.5943 

4.480 

0.920 

0.6943 

3.840 

0.851 

0.9820 

1.455 

0.524 
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Traverse  No.  6 


of  the  vertical 

traverse  from  the 

centre  =  3.S 

Distance  from 

Pressure  Diff. 

Velocity 

Upper  Wall,  in. 

Ah  cm  of  fluid 

u  ft/sec 

0.0080 

0.560 

0.3252 

0.0246 

0.590 

0.334 

0,0746 

0.940 

0.4215 

0.1246 

0.965 

0.427 

0.1746 

0.970 

0.4285 

0.2746 

0.955 

0.4250 

0.3746 

1,000 

0,435 

0.4746 

1.155 

0.467 

0.5666 

1.005 

0.436 

0.6666 

0.975 

0.4295 

0.8666 

0.395 

0,411 
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Table  9 


Flow  rate  of  water 

=21  USGPM 

Mean  Temperature  of  water 

=  75°F  (±2 

Reynolds  number  (based  on 

equivalent  diameter)  =  13650 

Traverse  No.  1 

Distance  of  the  vertical  traverse  from  the 

centre  =  0.000  . 

Distance  from 

Pressure  Diff, 

Velocity 

Upper  Wall.  in. 

Ah  cm  of  fluid 

u  ft/sec 

0.008 

1.485 

0,53 

0.0250 

1,685 

0.564 

0.0500 

2.967 

0.7475 

0.0750 

4,140 

0,884 

0,1750 

5,760 

1,043 

0,2750 

6,350 

1,095 

0,3750 

7.088 

1,157 

0.5000 

7.573 

1.196 

0.6250 

7,105 

1,159 

0.7250 

6,320 

1,092 

0,8250 

5,800 

1,047 

0.9250 

4.230 

0.895 
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Traverse  No.  2 


Distance  of  the  vertical  traverse  from  the  centre  =  1.000  in 


Distance  from 
upper  Wall,  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0080 

1.877 

0.595 

0.0250 

2.105 

0.631 

0.0500 

3.203 

0.777 

0.0750 

4.225 

0.894 

0.1750 

5,575 

1.027 

0.2750 

6.732 

1.131 

0.3750 

7.115 

1.159 

0.5000 

7,325 

1.177 

0.6250 

7.120 

1.159 

0.7250 

6.800 

1,132 

0.8250 

5.480 

1,019 

0.9250 

4.315 

0.902 

Traverse  No.  3 

Distance  of  the  vertical 

traverse  from  the 

centre  =  2 . < 

0.0080 

1.280 

0.492 

0,0250 

1.630 

0,555 

0.0500 

3,065 

0.751 

0.0750 

3,979 

0.867 

0.1750 

5.252 

0,996 

0.2750 

6,115 

1,073 

0.3750 

6 . 866 

1.139 

0.5000 

7.005 

1,15 

0.6250 

6.865 

1.139 

0.7250 

6.115 

1.073 

0.8250 

5.255 

0.996 

0.9250 

4.046 

0,874 
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Traverse  No.  4 


Distance  of  the  vertical  traverse  from  the  centre  =  3.000  in 


Distance  from 

Pressure  Diff. 

Velocity 

Upper  Wall,  in. 

Ah  cm  of  fluid 

u  ft/sec 

0.0080 

1.120 

0.460 

0.0250 

1.647 

0.558 

0.0500 

2.581 

0.699 

0.0750 

3.625 

0.827 

0.1750 

4.590 

0.931 

0.2750 

5.663 

1.035 

0.3750 

5.910 

1.058 

0.5000 

6.057 

1.070 

0.6250 

5.910 

1.053 

0.7250 

5.663 

1.035 

0.8250 

4.590 

0.931 

0.9250 

3.625 

0.827 

Traverse  No.  5 

Distance  of  the  vertical 

traverse  from  the 

centre  =  3.! 

0.0080 

0.845 

0.400 

0.0250 

1.232 

0.4825 

0.0500 

2 . 440 

0.679 

0.0750 

3.703 

0.836 

0.1750 

4.293 

0.901 

0.2750 

4.845 

0.956 

0.3750 

5.945 

1.062 

0.5000 

6.463 

1.106 

0.6250 

5.965 

1.062 

0.7250 

4.845 

0.956 

0.8250 

4.295 

0.901 

0.9250 

3.300 

0.846 
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Traverse  No.  6 

Distance  of  the  vertical  traverse  from  the  centre  =  3.750  in. 


Distance  from 
Upper  Wall,  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0080 

1.000 

0.434 

0.0250 

1.278 

0.491 

0.0500 

2.186 

0.643 

0.0750 

3.059  '■ 

0.759 

0.1750 

3.955 

0.364 

0.2750 

4,375 

0.909 

0.3750 

4.540 

0.925 

0,5000 

4.614 

0,933 

0,6250 

4.540 

0.925 

0.7250 

4.375 

0.909 

0.8250 

3.955 

0.864 

0.9250 

3.059 

0.759 
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Table  10 


Flow  rate  of  water  =  31.5  USGPM 
Mean  temperature  of  water  =  79°F(-2°f) 
Reynolds  number  (based  on  equivalent  diameter)  =  20,500 


Traverse  No.  1 


Distance  of  the  vertical  traverse  from  the  centre  =  0.000  in. 


Distance  from 
Upper  Walt  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0080 

2.022 

0.619 

0.0250 

4.146 

0.885 

0.0500 

7.415 

1.184 

0.0750 

9.805 

1.362 

0.1750 

12,326 

1.528 

0.2750 

13.915 

1.621 

0.3750 

15.182 

1.698 

0.5000 

15.754 

1.725 

0.6250 

15.132 

1,698 

0.7250 

13.912 

1.621 

0.8250 

12.348 

1.529 

0.9250 

9.738 

1.357 
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Traverse  No.  2 


Distance  of  the  vertical  traverse  from  the  centre  =  1.000  in 


Distance  from 

Pressure  Diff. 

Velocity 

upper  Wall  ,  in. 

/Vh  cm  of  fluid 

u  ft/sec 

0.0030 

2.862 

0,736 

0.0250 

4.2  30 

0.8945 

0.0500 

7.014 

1.150 

0.0750 

9.255 

1.323 

0.1750 

12.039 

1.510 

0.2750 

13.794 

1.615 

0.3750 

14.914 

1.630 

0.5000 

15.461 

1.713 

0.6250 

.14.946 

1.681 

0.7250 

13.800 

1.615 

0.8250 

12.043 

1.510 

0,9250 

9.456 

1.339 

Traverse  No.  3 

Distance  of  the  vertical 

traverse  from  the 

centre  =  2.( 

0.0080 

2.635 

0.7055 

0.0250 

4.255 

0,896 

0.0500 

7.289 

1,175 

0.0750 

8.863 

1.295 

0.1750 

10.988 

1.442 

0.2750 

12.611 

1.543 

0.3750 

14,-167 

1.638 

0.5000 

14.799 

1.673 

0.6250 

14.167 

1.638 

0.7250 

12.611 

1.543 

0.8250 

10.988 

1.442 

0.9250 

9.076 

1.309 
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Traverse  No.  4 


Distance  of  the  vertical  traverse  from  the  centre  =  3.000  in 


Distance  from 
Upper  Wall  ,  in. 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0080 

2.800 

0.728 

0.0250 

4.641 

0.936 

0.0500 

5.555 

1.023 

0.0750 

3.243 

1.162 

0.1750 

10.045 

1.380 

0.2750 

11.431 

1.472 

0.3750 

12.362 

1.529 

0.5000 

12.834 

1.559 

0.6250 

12.362 

1.529 

0.7250 

11.404 

1.470 

0.8250 

10.043 

1.380 

0.9250 

8.243 

1.162 

Traverse  No.  5 

Distance  of  the  vertical  traverse  from  the  centre  =  3.500  in 


0.0080 

1.595 

0.550 

0.0250 

3.446 

0.807 

0.0500 

5.464 

1.018 

0.0750 

6,353 

1.140 

0.1750 

8.938 

1,301 

0.2750 

10.217 

1.392 

0.3750 

11.160 

1.453 

0.5000 

11.806 

1.495 

0.6250 

11.160 

1.453 

0.7250 

10.217 

1.392 

0.8250 

8.966 

1.302 

0.9250 

7.006 

1.152 
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Traverse  No.  6 


Distance  of  the  vertical  traverse  from  the  centre  =  3.750  in. 


Distance  from 
upper  Wall 

Pressure  Diff. 

Ah  cm  of  fluid 

Velocity 
u  ft/sec 

0.0030 

1.811 

0.585 

0.0250 

3.507 

0.814 

0.0500 

5.656 

1.033 

0.0750 

7,589 

1.199 

0.1750 

9.263 

1.322 

0.2750 

9.596 

1.350 

0.3750 

10,484 

1,409 

0.5000 

10,693 

1,422 

0.6250 

10.484 

1.409 

0.7250 

9.596 

1.350 

0.8250 

9.285 

1.325 

0.9250 

7,589 

1.199 
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V  Basing  calculations  on  the  maximum  velocity  in  a  traverse, 
and  using  the  u”**  -  y'*"  correlation  of  Deissler,  the  friction 
velocity  u*  can  be  calculated.  Calculations  for  the  centre 
line  traverse  at  a  Reynolds  number  of  13,650  are  shown  here. 

Deissler 's  correlation 

u+  =  3.8  +  i —  In  y'*' 

0.360 


or 


=  3.8  + 


0.360 


In  ^ 


Solving  this  equation  for  u*  with 
u  =  1.196  ft/sec 

at  y  =  0.04167  ft. 

u*  =  0.06  ft/sec 

Hence  for  Traverse  No.  1  at  a  Reynolds  number  of  13,650 


Table  11 


u"^ 

13.6 

9.39 

27.2 

12.48 

40.7 

14,75 

95.2 

17.40 

149.7 

18.27 

204 

19.30 

272 

19.95 

Similar  calculations  for  Traverse  No.  2  at  the  same 
Reynolds  number  give  the  following  values  of  y***  and  u"*" 
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Table  12 


u 


+ 


13.6 

27.2 

40.7 

95.2 
149.7 
204 
272 


10.57 

12.95 

14.9 

17.12 

18.87 

19.33 

19.62 


The  values  given  in  Table  11  and  Table  12  have  been 
presented  in  Figure  13  and  Figure  14  respectively. 
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FIGURE  2 


THE  TRAVERSING  MECHANISM 
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FIGURE  3  u+-y+  Plot  of  the  Data  on  Tube 


0.961  Umax, 
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FIGURE  5  Velocity  Distribution  0.000  in.  from  centre  at 
Reynolds  Number  of  13650 
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Velocify,  ft  / sec. 

FIGURE  6  Velocity  Distribution  1.000  in.  from  centre  at 
Reynolds  Number  of  13650 
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Velocity,  ft/  sec  . 

FIGURE  7  Velocity  Distribution  2.000  in.  from  centre  at 
Reynolds  Number  of  13650 
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Velocily,  f+/sec. 

FIGURE  8  Velocity  Distribution  3.000  in.  from  centre  at 
Reynolds  Number  at  13650 
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Velociiy,  f+/sec, 

FIGURE  9  Velocity  Distribution  3.500  in.  from  centre  at 
Reynolds  Number  at  13650 
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Velocity,  ft  / sec. 

FIGURE  10  Velocity  Distribution  3.750  m.  from  centre  at 
Reynolds  Number  of  13650 


84 


saqouT  'aqeid  uioq.qoq  mojLj  qqbxaq  xauueqo 


0.5  I  1.5  2  2.5  3  3.5 

Channel  Width  from  Centre,  inches 

FI  CURE  II  Cons  +  an+  Veloci  +  y  Lines  at  a  Reynolds  Number  of  13650 
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FIGURE  13  Plot  of  Centre  Line  Traverses  Data, 
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FIGURE  14  u  -  y  +  Plot  of  Data  on  Traverses  I  in,  from  the  centre 
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FIGURE  15  Approximate  Shear  Stress  at  the  Wall. 
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NOMENCLATURE 


AA 

b 

c 

d 

E 

f 

H 

L 

AL 

P 

Re 


S 

t 

u 

U 


u* 


area  between  two  velocity  gradient  lines 
half  the  height  of  the  channel 
velocity  of  light 
diameter  of  tube 

e.m. f . 

friction  factor 
gravitational  constant 
magnetic  field  strength 
length  of  tube 

distance  between  velocity  gradient  lines  at  the  wall 

pressure 

Reynolds  number, 

distance  between  electrodes 
time 

point  velocity 

average  velocity,  also  used  for  point  velocity 
in  universal  velocity  distribution  equations 
dimensionless  velocity  parameter,  u/U* 

1/2 

friction  velocity,  (Tq  9q/p) 


V 


velocity 
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w 

U  '  ,  V  '  ,  w' 


X,  y,  z 


velocity 

fluctuating  velocity  components 
bar  indicates  time-average  values 
cartesian  co-ordinates 


T 

P 

-t) 


shear  stress 
density 

dynamic  viscosity 
kinematic  viscosity 
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